Draft version September 28, 2010 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



O 

D 
C/^ 

l> 
(N 

o 

(N 
> 

ON 

O 
O 



ACCELERATION OF ULTRA-HIGH ENERGY COSMIC RAYS IN THE COLLIDING SHELLS OF BLAZARS 
AND GRBS: CONSTRAINTS FROM THE FERMI GAMMA RAY SPACE TELESCOPE 

Charles D. Dermer^ & Soebur Razzaque^'^ 

^ Space Science Division, Code 7653, Naval Research Laboratory, 
Washington, DC 20375-5352, USA and 
^NRL/NRC Research Associate 

(Received April 23, 2010) 
Draft version September 28, 2010 

ABSTRACT 

Fermi Gamma ray Space Telescope measurements of spectra, variability time scale, and maximum 
photon energy give lower limits to the apparent jet powers and, through 77 opacity arguments, the bulk 
Lorentz factors of relativistic jets. The maximum cosmic-ray particle energy is limited by these two 
quantities in Fermi acceleration scenarios. Recent data are used to constrain the maximum energies 
of cosmic-ray protons and Fe nuclei accelerated in colliding shells of GRBs and blazars. The Fermi 
results indicate that Fe rather than protons are more likely to be accelerated to ultra-high energies 
in AGNs, whereas powerful GRBs can accelerate both protons and Fe to > lO^'^ eV. Emissivity of 
nonthermal radiation from radio galaxies and blazars is estimated from the First Fermi AGN Catalog, 
and shown to favor BL Lac objects and FRl radio galaxies over flat spectrum radio quasars, FR2 
radio galaxies, and long-duration GRBs as the sources of UHECRs. 

Subject headings: cosmic rays - galaxies: active - galaxies: jets - gamma rays: galaxies - radiation 
mechanisms: nonthermal - shock waves 



1. INTRODUCTION 



Mpc in order to reach us without los in g significant 



'Hillas' (1984") pointed out an essential requirement for 
acceleration of ultra- high-energy cosmic rays (UHECRs), 
namely that the particle Larmor radius 7'l = E/QB must 
be smaller than the size scale of the acceleration region. 
Here E is the particle energy, Q = Ze is its charge, and 
B is the magnetic field in the acceleration zone. This 
limitation applies to Fermi acceleration scenarios where 
a particle gains energy while diffusing through a mag- 
netized region. Additional limitations due, for exam- 
ple, to radiative losses or available time, further restrict 
the maximum energies and therefore the allowed sites of 
UHECR acceleration. 

Two plausible classes of astrophysical accelerators 
of extr agalactic UHECRs are acti ve galactic nuclei 
(AGN ) (jMannheim fc BiermaimI 1198 9: Bcrezinsk v et al. 
2002) and gamm a-ray bursts (GRBs) (jWaxman 



19951 iVietr i l il995l : [Milgrom fc Usov 1995) (see also 
Rachen fc M cszaros 1991 [Halzen fc Hoopeil 120021: 



Dermer fc Menon.2009i) . though other types of sources, 
including young, hi ghly magnetized neutron stars 
(Ghisellini et al.l 120081 ) and structure formation shocks 



jlnoue 20dl" remain viable . The announcement by the 
lAuger collaboration! (|2007[ ) of anisotropy in the arrival 
directions of cosmic rays with energies E >6 x 10^^ eV, 
even given the reduced correlation in the latest data 
from the Pierre Auger Observatory ([Abraham et al.l 
I2009D . is compatible with the production of UHECRs 
in many source classes, including GRBs and blazars. 
Because of the GZK effect involving photohadronic 
interactions of protons or ions with CMB radiatio n 
(IGreisenlll966t IZatsepin fc Kuz'minl[T96l ISteckei<[T968l) . 
higher-energy cosmic rays with E > 10^° eV must be 
produced by sources located within distances d < 100 

[charles . dermer @nrl . navy, mil] 



energy (e.g., INagano fc Watsoiil I2000t IHarari et al.l 
2006). 

The most powerful AGNs and long-duration GRBs are 
found far outside the GZK radius, at redshifts z > 0.1. It 
is therefore of interest to re-examine Fermi acceleration 
requirements to determine if there are AGN and GRB 
sources with appropriate properties within the GZK ra- 
dius. Here we make a detailed examination to justify 
a simple derivation of maximum particle energy relat- 
ing apparent source power and bulk Lorentz factor F in 
the framework of Fermi acceleration in colliding shells. 
(Note that these arguments do not apply to non-Fermi 
type mechanisms, for example, electrodynamic accelera- 
tion in pulsar magnetospheres.) The derived limits are 
compared with values implied by Fermi data, yielding 
constraints on UHECR acceleration in these sources. We 
then use the First Fer mi Large Area Tele scope (LAT) 
AGN Catalog (ILAC) (|Abdo et al.ll2010aD to estimate 
the nonthermal emissivity of AGNs. 

We find that the lower luminosity BL Lac objects and 
FRl radio galaxies are more likely to be the sources 
of UHECRs than the rare, powerful flat spectrum ra- 
dio quasars (FSRQs) and FR2 radio galaxies, and are 
more likely to accelerate Fe than protons to ultra-high 
energies. GRBs, on the other hand, can accelerate both 
protons and Fe nuclei to ultra-high energies, but are rare 
within the GZK volume. 

2. MAXIMUM PARTICLE ENERGY IN COLLIDING SHELLS 

The total comoving energy density u', including rest- 
mass and magnetic-field energy density, of a cold, 
isotropic relativistic wind with total wind power L and 
outflow Lorentz factor F = 1/^1-/32 at radius r from 
the source is u' — L/ (47rr^/3r^c). Primes here and below 
refer to quantities measured in the proper (comoving) 
frame of the radiating fluid. If a fraction of the to- 
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tal energy density is in the form of magnetic-field energy 
density u'g, = i?'^/87r, where B' is the magnetic field in 
the fluid frame, then rB'V = ^JiesL / j3c, implying maxi- 
mum particle energies E^ax — pVQB'r' = (3QB'r (since 
the comoving size scale r' = r/T). Thus 



eV 



(1) 



where the nonthermal 7-ray luminosity — 
10''^L46 erg s~^ (e.g.. lWaxmanll200l:lFarrar fc Gruzinovl 
dOOl), and we write = CgL, where eg is the fraction of 
jet energy in electrons that is assumed to dominate the 
radiative 7-ray output. In general, Lj < L, and <^ L 
for radiatively inefficient flows. Besides giving the min- 
imum apparent isotropic jet power, 7-ray observations 
give minimum values of F from 77 opacity arguments, 
allowing us to identify whether a given source is a plau- 
sible site for UHECR acceleration. 

The estimate in eq. ((!]) does not, however, explain 
how a cold magnetohydrodynamic wind can transform 
directed kinetic energy to relativistic particles, which 
requires consideration of a specific model. Within the 
colliding shell framework (Recs & Mcszaros 1994; Piran 
I1999D . which is often invoked to explain the formation 
of spectra in GRBs and blazars, we can assess the con- 
ditions under which eq. ([T|) is valid. In this model, a 
central black hole is assumed to eject a slower shell a 
with coasting Lorentz factor Fq — F^ during explosion 
frame times < < At^a, followed by a faster shell b 
with Fo = F{, > Fa ejected at times i*d < t* < i* -I- Ai^^t, 
where t^d is the stationary-frame delay time between the 
start of the ejections of shell a and b. The shell ener- 
gies £a(b) are related to their luminosities L^^f,) through 

£a{b) = La(h)A't*a(b)- 

The shells are assumed to collide after they reach their 
coasting phase. Neglecting shell spreading (which can 
be included by renormalizing the shell durations), and 
assuming that the event takes place sufhciently quickly 
so that we can approximate the shell density as constant 
during the duration of the collision, then simple kinemat- 
ics shows that the collision radius Vcoii and collision time 
t*,coii are given, in the limit F^ ^ 1, by 



'^coll — ct^ 



Pt 



t b 



(2) 

The proper frame particle density in the shells is 
"■l(fc) = La(b)/'^T^^'i(pf'^rnpC' at radius r. Letting F(> 1) 
denote the Lorentz factor of the shocked fluid, then the 
Lorentz factors of the forward / and reverse r shocks 
as measured in the shocked fluid frames are F 

FF 



/('■) 



aX6)(l-/3a(&)/3) ^ W/^a{b)+ya{b)m (|Sari fc Piran| 



1995!). From the equality of kinetic energy densities at 
the contact discontinuity, we have 



f2 _ r 



FH Fr F^ La La 



Four asymptotic regimes can be identified, depending 
on whether the forward shock is relativistic (RFS) or 
nonrelativistic (NFS), and the reverse shock is relativistic 
(RRS) or nonrelativistic (NRS): 



u = n' In' 

b a 




4. NRS 
NFS 



r/r 



Fig. 1. — Different regimes in eoUiding shell interaetions. 

1. RRS (F^ > 1) and RFS (F/ > 1). This holds 
when Ffc > F > Fq , implying 

,1/4 



F = ^1/4, /p r, Ff = , and F^ = — ttt 



2ui/4 

when pr <C mm{^/u, 1/y/u). 
2. NRS (F^ - 1 « and RFS, implying 



(4) 



F - Fb , F/ = , and /3. = (5) 



when 1/v2m <C pr ^ 1/2 or ?i ^ 1. For this case. 
Lb > La. 



3. RRS and NFS (F/ - 1 w /3j/2), implying 

^-^'^'^^ = ^/I'""^^'--2^ 
when ^/u/2 < pr < 1/2 or u < 1. 



(6) 



4. NRS and NFS, that is, < 1 and I3f < 1, which 
takes place when Fa = Ff, = F. Because pr — 1, 
the duration of the interaction for this case may be 
sufficiently long to violate the assumption of con- 
stant shell density during the collision, though we 
formally treat it here. 

Fig. 1 illustrates the various regimes in terms of the 
parameters u and pr for which we derive maximum par- 
ticle energies. The magnetic field B ' in the shocke d 
fluid is defined in the usual way re.g.. lSari etHI [T998h . 
through a magnetic-field parameter es, so that B^j^) = 

^327rTOpc2n^^^^eB/(r)(F^j,^^ - Tf(r))- The duration of 

the flare from the forward shock is determined by the 
time At'a required for the forward shock to pass through 
shell a. The width of shell a in the frame of the shocked 
fluid is Aa/F/ due to length contraction, where the 
proper frame width of shell a is Aa = FAa(r). If /3/c is 
the speed of the FS, then /?/ ^ 4/3//3 when F/ - K 1, 
and /3/ ^ 1, f/ = \/2F/, when F/ > 1. The duration 
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of the event is also limited by the possibility that shell 
b dissipates its entire energy before the forward shock 
passes through shell a. Using similar reasoning for the 
reverse shock, we obtain the comoving timescales for the 
interactions at the forward and reverse shocks: 



At', 



FS{RS) 



^ ^a{b)^a{b){r) 

/3/wr/(r)c 



-6(a) 



(7) 



Following the passage of the reverse and forward shocks 
through the shells, adiabatic expansion quickly ends sub- 
sequent acceleration and emission (Dermer ,2008i ). 

The maximum energy of particles accelerated at the 
forward and reverse shocks is given by E^axjir) — 

ZeTB'j:f^^^cPf(^r)^t'ps{RS)- general expression can be 
written as 



'max,f{r) 



Ze 



2eBf{r)La(b) f r .2 



1 



-fir) 



(1 -p2)min[l. 



r(r/M-i)J 



(8) 



The general expression for the maximum radiative effi- 
ciency giving the internal energy dissipated in the for- 
ward and reverse shocks can be written as 



Ef(r) = min [Sa(b) 



r(r/(,.) - 1) 



a{b) 



(9) 



For the different cases, we obtain the following results: 

2.1. Acceleration at the Forward Shock 
1. RRS and RFS. 



Ejnax^f — ^ 



Ze _ j2tBfLb 



min[l. 



£f = min[£b. 



{tn^d At^a) I At^a 
£a Lb, 




(10) 



2. NRS and RFS. 



r 



Ze hcBfLb jL^£b2miTi[l,£a/2pl£b] 



Lb £a (t*d Attfa) I Ai^^ 

mm[£b,£a/2pl] . (11) 



3. RRS and NFS. 



Ze 2eBfLb mm[l, 4{£b/£a){La/ Lb] 



£ 



f 



C {tf.d — At^a) / A.t^:a 

UYm[£aA£b{La/Lb)], (12) 



4. NRS and NFS (also written for acceleration at the 
reverse shock). 



max ,f {r) 



Ze _ j2eBfLa{b) 



V2 



{t^d ~ At^a)/At^a 



^f(r) = min[/3/(r)^-a(b)/2,'?b(a)] • 

2.2. Acceleration at the Reverse Shock 
1. RRS and RFS. 



(13) 



Er, 



Ze 2tBrLb 2min[l,£6/2f,] 



r V C [t^d - At^a)/^t*a 

£r^mm[£a.£bl2]. (14) 



2. NRS and RFS. 



Ze [2^^b [l: 2min[l,^] 



r 



Lb {t*d — At^a) / A.ti,a ' 

£r = mii\[£a,£b/{2upY)] , (15) 
noting that up^ ^ 1 for this case. 
3. RRS and NFS. 



Ze / 2eBrLb min[l, 2£a/£b\ 

~y ~C {t^d ~ At^a)IAt^a ' 

£r = min[£a,'^^6/2] . (16) 



En 



4. NRS and NFS, given by eq. 

3. LIMITATIONS ON UHECR ACCELERATION 

The maximum particle energy for the various cases is 
always proportional to the umbrella function, eq. ([T]), 
derived from elementary principles, but multiplied by a 
coefhcient < 0(1). The ability of a shell collision to ac- 
celerate particles to the highest energies is conditioned 
on very specific behaviors of the shells, namely that the 
second shell is much faster than the first (pr ^1): and 
that the time t^,d between shell ejections is a small factor 
larger than the duration At^a of the event forming shell 
a (as expressed by the term [t^d ~ At^a) / Ati,a in the 
denominators of eqs. - p6)) ). The most favorable 
regime for particle acceleration to the highest energies 
occurs for the case of a RRS and RFS when the energies 
and luminosities of the two shells are about equal. This 
also gives the highest radiative efficiencies. The main 
requirement is a larg e contrast b etween t he Lorentz fac- 
tors o f the two shells ()Beloborodo v 2000; Kumar fc PiranI 
[2000). 

The highest radiative efficiency coincides with approx- 
imately equal energies and luminosities for the cases of 
a RRS and RFS, and a RRS and NFS. In the case of a 
NRS and RFS, where Lb >• La is required for validity of 
this asymptote, a much larger energy in shell h than shell 
a is required for maximum radiative efficiency at the re- 
verse shock, as shown by eq. (|15p . Energy dissipation 
in this case would, however, more likely be dominated 
by the forward shock. Kinematic limitations ensure that 
the radiative efficiency is poor for dissipation at either 
the forward or reverse shocks for the case of a NRS and 
NFS, eq. ([T3| . depending on the precise energies in each 
of the shells. 

Supposing that the engines of GRBs and blazars or, for 
that matter, microquasars, eject shells with such prop- 
erties (which is necessary in the case of GRBs to explain 
their high 7- ray radiati ve efficiency in an internal shock 
scenario; cf. Iloka et aLii2006, ). then we can construct a 
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diagram illustrating the viability of various sources to ac- 
celerate UHECRs. In Fig. 2 we plot apparent luminosity 
as a function of Lorentz factor for the acceleration of lO^*^ 
eV protons (heavy solid curve) and Fe nuclei (heavy dot- 
dashed curve), from eq. Note that we plot this equa- 
tion to nonrelativistic values of F, which is outside the 
regime where we have considered its validity in a colliding 
shell scenario. Sources above these curves can in princi- 
ple accelerate UHECRs. The acceleration rate in Fermi 
scenarios is governed by the Larmor timescale, so that 
the acceleration timescale t'^^c — (f'^'/c = (pE' /ZeB'c, 
and (/> ^ 1. Incidentally, the requirement that t'^^c is 
smaller than the available time t'^^^ = Tty / {1 + z) , where 
ty is the measured variability timescale, essentially recov- 
ers eq. (dl when 0=1 and r = T"^ ct^ / + z) . This shows 
that eq. ([T]) is a restatement of the Hillas condition by 
relating B to L and tl to ty. 

Maximum particle energy is also limited by the require- 
ment that (i) t'^ec is short er than the synchrotr on energy- 
loss timescale t'^y^ (e.g.. iGuilbert et aLlll983() . We can 
also consider the (less restrictive) condition (ii) t'^ 



syn 



> 



t'g^ya, so that particles do not significantly cool during 
the available time. Writing the comoving magnetic field 
u'b = eBL^/{ee^Tir^T^c) and E20 = E/lO'^^eV, the for- 
mer constraint becomes 



rp a2 



.^3/2 



Q'Kecty 



where Amp is the particle mass, implying 



(17) 



(18) 

For the second case, 
E 



— = 37r(-j [—) — -r-r TT-T ' (19) 



a^L^ 1 -I- z 



implying 



^-^^"'^^ £;.o(i + z) ^''^ 

The restrictions implied by eqs. and ([20)1 are 

shown by the dashed and dotted lines, respectively, in 
Fig. 2, for parameters characteristic of UHECR pro- 
ton acceleration to VS^^ eV in blazars and GRBs. Here 
^e/^B = ^i'i> = 10, and ty = l{f s and 10 ms, and F = 10 
and 10^, for blazars and GRBs, respectively. We also 
plot data for various sources observed with Fermi and 
ground-based 7-ray telescopes. In all cases except Cen- 
taurus A, T^in is derived from 77 opacity arguments, 
with the apparent 7-ray luminosity giving the minimum 
source luminosity. The inference of Tmin from 77 opac- 
ity arguments is model dependent, with the determina- 
tion of T^in dependent on assumptions about target- 
photon anisotropy, relationship between variability time 
and emission region size scale, photon escape probabil- 
ity, and the dyna mic state of the emitting plasma (e.g., 
lAckermann et al.ll2010l:lGranot et ai]|2008D . Even so, the 
strong F-dependence of comoving photon energy density 
u' (X makes it unlikely that the actual value of Fmm 



differ by more than a factor of w 2 from the value derived 

through simple 77 arguments. 

The long-duration GRB 080916C (lAbdo et al.ll2009aD 
and the short-duration GRB 0905 lOA (jAckermann et aLl 



IMof) have Tmin ~ 10^ and w 10^^ gj.g g-i Qp-g 
090902B, with Tmin ~ 10^ and L-y w 10^^ erg s~^ be- 
tween 6 and 13 s after the trigger time (jAbdo et al.l 
l2009bD . would cluster in the same regime. For 3G 454.3, 
Trmn « 8 and L.^ « 5 X 10** erg s'^ fAbdo et al.l[2009G ). 



In the case of NGC 1275, 



10^2 erg s ^ and the 



Doppler factor (and therefore F) is > 2 (jAbdo et al.l 
[2009d). For PKS 2155-3 4, a B L Lac object, we use 
the results of Finkc ct al. (2008") for the giant flares of 
2006 July (Aharonian ct al. 2007), which employs a syn- 
chrotron self-Gompton model with 77 absorption and 
various EBL models to derive F « 100. Note that the 
absolute jet powers derived there can be much less than 
the apparent jet power that enters into the acceleration 
constraint defined by HESS measurements of its appar- 
ent isotropic 7-ray luminosity w 3 x 10^^ erg s~^. 

Finally, we consider the case of the FRl radio galaxy 
Cen A, the only one of the sources shown in Fig. 1 that is 
within the GZK radius. It is of special interest, of course, 
because of the clustering of t he arrival directions of sev- 
eral UHECRs towards Ccn A ( Auger collabor ationll2007l: 
[Moskalcnko ct__aj, 2009; Abraham ct al. 2009), and early 
speculatio ns that it could b e a dominant source of UHE- 
CRs (Piran &: Farrad[2001h . Because its jet is pointed 
away from our line of sight, the jet luminosit y and F fac- 
tor o f Cen A can only be indirectly inferred ( Kraft et al.l 
l2002f ). One way is to assume that the energy of the ra- 
dio lobes are powered by the jets, and use synchrotron 
theory and lobe dynamics to infer total energy and life- 
time. Values between 10^^ _ iQ^^ gj-g g-i are inferred 
(jHardcastle et al.ll2009f ). with jet beaming and episodes 
of intense outbursts arguably capable of allowing the jet 
to reach apparent powers sufficient to accelerate UHECR 
protons (D ermer et a l. 2009). 

Deceleration of relativistic shells by the surround- 
ing medium generates a relativistic external for- 
ward shock and a relativistic/non- relativistic reverse 
shock. The time scale for deceleration depends on 
the apparent isotropic kinetic energy of the merged 
shells -Ek.iso, the bulk Lorentz factor Fq and the 
density of the surrounding medium n, given by 
tdec = (1 + z)(3^k,iso/[32^mpc5nFg])i/3 « 1.9(1 -f 
z)n~^/3£;i/3p-s/3 rpj^^ subsequ ent evolution of the 
blast wave is described by the iBlandford &: McKei 
(|1976f ) self-similar solutions. Acceleration of cosmic rays 
to maximum energies in the forward shock takes place 
during the deceleration time, and similar to eq. ([T]) we 
can write from i^cc = iacc 



Ze F^/^ 

En... = — -— ^(97r£;2iso4"™pc2)V6 



1.4 X 10 



21 



(1 + ^) 



n^/'e]i'E}fTr eV. (21) 



The corresponding constraint on the apparent isotropic 
kinetic energy and bulk Lorentz factor to accelerate par- 
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tides to 10^° eV is 

i?k.iso « 3.4 X 1054^!ii±i erg. (22) 

This constraint is satisfied by most long-duration GRBs 
to accelerate both protons and Fe nuclei, and by FSRQ 
blazars (depending in detail on the energy of the blazar 
flare and density of decelerating medium) to accelerate 
Fe nuclei. 

4. DISCUSSION 

Fig. 2 shows that the short- and long-duration GRBs 
for which Fermi observations give both and Tmin 
easily satisfy the luminosity requirements to accelerate 
UHECR protons or ions. After considering the specific 
parameter values that enter into eqs. (IT51) and (PO)) . one 
finds that the additional constraints imposed by the syn- 
chrotron cooling rate are not severe either for blazars 
or GRBs. One difficulty for arguing that GRBs are the 
sources of UHECRs is their rarity within the GZK ra- 
dius. Only if the intergalactic magnetic field is suffi- 
ciently strong (^ nG with Mpc scales for magnetic-field 
reversals) to disperse the arrival time of the UHECRs, 
but not so strong to erase their inhomogeneous arrival di- 
rections, ca n long and short GRBs be plausible UHECR 
candidates ([Razzaque et al.l I2009D . More complicated 
magnetic field g eometries can also relieve this problem 
(|Kashti fc WaxmanI 12001) . A further difficulty accom- 
panying the large F values and correspondingly dilute 
comoving photon energy densities implied by the Fermi 
results on GRBs is that photohadronic processes are sup- 
pressed. Intermediate neutron production with the es- 
cape of ultra-high energy neutrons that subsequently de- 
cay to form UHECRs was prop osed as a principal mech- 
anism (jAtovan fc Dermeill2003D to circumvent the prob- 
lem of the escape of UHECR ions. Such escape is prob- 
lematic in a bursting source because the strong flux of 
ions will generate a shock that causes the ions to lose en- 
ergy as they leave the GRB. (Murase fc Beacom 2010). 
Future Fermi observations will reveal whether there is a 
large population of low F-factor GRBs, or i f another class 
of GRBs, such as low-luminosity GRBs (|Murase et al.l 
[20061 IWang eFaHIMTt ILiang et al.l[2007l) . can make the 
UHECRs. 

From the ILAC catalog (lAbdo et al.ll2010af) . we can 
make a diagram. Fig. 3, of the volume-averaged nonther- 
mal 7-ray luminosity density (or emissivity). Here we 
use the time-averaged 100 MeV - 100 GeV luminosity 
measured over eleven months, and divide by the volume 
I'd associated with the proper distance d of the indi- 
vidual sources to make a cumulative emissivity for differ- 
ent classes of 7-ray galaxies. The volume-averaged emis- 
sivity, unlike the source density, is independent of the 
beaming factor. The cumulative emissivities are shown 
separately for BL Lac objects, FSRQs, misaligned AGNs, 
and non-AGN star-forming and starburst gala xies, in- 
cluding M82 and NGC 253 (lAbdo et al.ll2010bl) . as weU 
as NGC 4945 reported in the ILAC. NGC 4945 is classi- 
fied here as a starburst, though it also contains a Seyfert 
nucleus. The misaligned AGNs consist of 11 sources, 
including seven FRl radio galaxies and four FR2 ra- 
dio sources. The FRl galaxies are Cen A, M87, NGC 
1275, NGC 6251, NGC 1218 (3C 78), and PKS 0625-35 



(|Abdo et al.ll2010af) . The FR2 objects consist of two ra- 
dio galaxies, 3C 111 and PKS 0943-76, and two steep 
spectrum radio quasars, 3C 207 and 3C 380. 

For comparison with the cumulative emissivity, the 
fiducial luminosity-density value of « lO''^ erg Mpc"^ 
yr^^ that is needed for clas ses of sources to energize 
UHECRs against GZK losses (|Waxman fc Bah call 199^ 
is shown. What is obvious from Fig. 3 is that FSRQs 
do not have sufficient emissivity to power the UHE- 
CRs under the assumption that the 7-ray luminosity 
is a good measure of the UHECR power. A much 
larger energy release in UHECRs than 7 rays is pos- 
sible, but even so, FSRQs are absent within the GZK 
radius, and FR2 radio galaxies, which are the putative 
parent p opulation o f FSRQs under the unification hy- 
pothesis ()Urrv fc Pa dovani 1995 ), are onl y found at dis- 
tances > 100 Mpc (Moskal enko et al.l !2009'). Pictor A is 
the closest FR2 radio galaxy at z = 0.035, and Cygnus 
A is at z = 0.056; neither has yet been reported as Fermi 
LAT sources. The redshifts of the detected FR2 radio 
galaxies 3C 111 and PKS 0943-76 are 0.049, and 0.27, 
respectively. 

The comparison is more favorable for BL Lac ob- 
jects which, as indicated by Fig. 3, have the necessary 
nonthermal power to energize UHECRs. The ones de- 
tected at GeV energies are still outside the GZK radius, 
though the famous TeV (and GeV) blazars Mrk 421 and 
Mrk 501 reside, at » 130 Mpc, just outside it (see also 
iTakami fc S ato 2009. for a discussion of the space density 
of putative UHECR sources) . The unification hypothe- 
sis would then suggest that many FRl radio galaxies are 
found at closer distances, including misaligned galaxies 
detected at 7-ray energies. Indeed FRl galaxies detected 
at GeV energies are found within the GZK radius, as 
shown by the misaligned AGNs in Fig. 3. Misalignment 
means that the 7-ray luminosity when viewed directly 
along the jet is probably much larger than the luminos- 
ity measured with Fermi. NGC 1275 has apparent 7-ray 
luminosity of ss 10^"* erg s~^ and, at a dis tance of w 75 
Mpc, faUs within the GZK radius (|Abdo e t al. 2009d). It 
is variable at 7-ray energies, indicating that much of its 
7 ray flux is probably associated with a jet. Moreover, 
it is a compact symmetric object, with transient out- 
bursts with durations of « 10 - 10^ yrs during which 
conditions are more favor able for UHECR acceleration 
(jHoriuchi fc Takamil 120091 ). The variability is not short 
enough to give Tmm from 77 arguments, though mod- 
eling results and observations of apparent superluminal 
motion suggest mildly relativistic Lorentz factors. By 
comparison, Cen A has relatively small apparent 7-ray 
lobe and core luminosities, each amounting to ~ 10^^ erg 
s^^, but its emissivity is large due to its proximity. 

The star-forming and starburst galaxies are abundant 
within the GZK radius, and have substantial 7-ray emis- 
sivity, which is more than adequate to account for the 
power needed to accelerate UHECRs. Where this source 
class falters, however, is in the low individual GeV - 
TeV 7-ray luminosities, representing « 3 x 10'^® erg s~^ 
for the Milky Way, and ~ 10'^° e rg s^^ for the lumi- 
nous star bursts (jAbdo et al.ll20ldbl ). Fermi acceleration 
of UHECRs with such low powers is, as seen from Fig. 
[21 not feasible. Furthermore, the lack of reported detec- 
tion of 7-ray emission from clusters of galaxies weakens 
the case for UHECR acceleration in structure-formation 
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shocks. 

Long-duration GRBs have apparent 7-ray luminosities 
far greater than needed to accelerate cosmic-ray protons 
or ions to E > 10^° eV, as shown in Fig. 2. Their time- 
averaged photon luminosity density is, however, insuf- 
ficient to power UHECRs within the GZK radius un- 
less the typical baryon loading in GRBs, representing 
the ratio of energy in cosmic rays to that radiated as 
photons, is ^ 1. Estimates for the local {z <^ 1) lu- 
minosity density, based on the luminosity function and 
local event rate density of long-duration GRBs, range 
from « 6x 10'*2(At/20 s) erg Mpc'^ vr-^ IScl5^ (pOOll 
though without using GRB redshift information), to ~ 
2x 10 '^'^(At/10 s) erg Mpc^^ yj.-i (|Wanderman fc PiranI 
|2010( ). where At is the mean duration of long GRBs in 
the explosion frame. A local luminosity density of (5 - 
8) X 10'^' ^(Af/10 s) erg Mpc" '^ yr~Ms derived in the treat- 
m ent of | Guetta et al.l (|2005l ) and the luminosity function 
of iLiang et al.l (|2007D implies a local luminosity density 
of« 2x 10'*3(Ai/10 s) e rg Mpc~^ yr~\ Based on a phys- 
ical model of GRB jets, iLe fc Permed (|200 f) calculate a 
long duration GRB density of (3 - 4) x erg Mpc~^ 
yr"-'^, assuming At = 10 s. 

By comparison with the nonthermal luminosity den- 
sity of long-duration GRBs, that of FRl radio galax- 
ies and BL Lac objects within the GZK radius is at 
least 1-2 orders of magnitude larger (Fig. 3). If the 
comparison is with the nonthermal emission emitted in 
the GeV/LAT range rather than at MeV energies, which 
could be thermal or photospheric emission, then the re- 
quired baryon loading must be an order of magnitude 
larger ()Eichler et al.ll20T9 ). The local ph oton luminos- 
ity densities of the short hard GRBs fe.g., |Guettall2006D 
or low luminosity GRBs (iWang et al.ll2007t ILiang et al.l 
[2OO7I: IMurase et all [200I " can also be comparable to 
the emissivity from long-duration GRBs, though with a 
larger local space density and smaller energy release per 
event. 



5. CONCLUSIONS 

Fermi observations shown in Fig. 3 indicate that FRl 
radio galaxies and misaligned BL Lac objects located 
within the GZK radius have sufficient emissivity to power 
the UHECRs. With typical Lorentz factors w 2 - 10, and 
apparent jet powers ss 10^** - 10^^ erg s~^ (which could 
exceed lO**^ erg s~^ and large Lorentz factors during flar- 
ing episodes). Fig. 2 shows that acceleration of Fe nuclei 
in FRl radio galaxies is possible in colliding shells made 
in the jets of these galaxies. Given the favorable circum- 
stances needed for colliding shells to accelerate UHECRs, 
including large Lorentz factor contrast and short times 
between shell ejections, the acceleration of protons is less 
likely. The L — T constraint is unfavorable for UHECR 
acceleration at sites with low apparent luminosity, such 
as starburst galaxies or the lobes of radio galaxies. 

Long-duration GRBs have sufficient power to acceler- 
ate cosmic rays to ultra-high energies, but their local 
photon luminosity density in photons, ~ 10'*'^ - lO**"* erg 
Mpc^'^ yr"^, implies comparable or large baryon loading 
in most models for UHECR acceleration. When com- 
pared with the clearly nonthermal Fermi LAT flux, the 
re quired baryon - loadin g becomes significant, as shown 
bv lEichler et al.l ()2010t ). The local nonthermal luminos- 
ity density of FRl radio galaxies and BL Lac objects by 
far dominates that of GRBs, especially when compared 
only with the LAT fluxes from GRBs and blazars. This 
circumstance favors UHECR acceleration by the super- 
massive black-hole engines in radio galaxies and blazars, 
provided that UHECRs are predominantly Fe ions. 
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Fig. 2. — Sources with jet Lorentz factor F = \//3^F^ + 1 must have jet power L exceeding the heavy solid and dot-dashed curves to 
accelerate protons and Fe, respectively, to i? = lO^'^ eV, from eq. JTJ. Upper limits to L as a function of F for acceleration of UHECR 
protons to lO'^" eV in blazars and GRBs are given by the dashed lines due to competition between synchrotron losses and acceleration, and 
by the dotted lines when comparing synchrotron losses and available time. Here we use variability time t„ = 10'' s and ^ ■C 1 for blazars, 
and t„ = 10 m s and z = \ for GRBs, as labeled, with = 10 in both cases. Scalings for different values of ti,, F, 2, and </> are given by eqs. 
(Hg and l(20) l. 
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Fig. 3. — Nonthermal luminosity density of different classes of 7-ray galaxies detected with Fermi. Here we show cumulative emissivities 
for FSRQ, BL Lacs, FRI and FR2 radio galaxies, and star-forming galaxies. The band between 100 and 200 Mpc labeled "GZK" 
represents the outer perimeter from which sources of UHECRs with energies > 10^" eV can originate. 



